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Abstract. In this paper we consider a model for a nucleon interacting with 
the uj and a mesons in the atomic nucleus. The model is relativistic, but 
we study it in the nuclear physics nonrelativistic limit, which is of a very 
different nature from the one of the atomic physics. Ground states with a 
given angular momentum are shown to exist for a large class of values for the 
coupling constants and the mesons' masses. Moreover, we show that, for a 
good choice of parameters, the very striking shapes of mesonic densities inside 
and outside the nucleus are well described by the solutions of our model. 



1. Introduction 

This article is concerned with the existence of ground states for a stationary 
relativistic mean-field model for atomic nuclei in the nuclear physics nonrelativistic 
limit. Moreover, a discussion about some interesting qualitative properties of the 
solutions of the model considered is provided. 

Classically, at low energies, nuclear structure is described by quantum mechanical 
many-body problems of fermions interacting by a nonrelativistic interaction (see for 
instance [7]), and as explained in [BJ, this interaction is understood to have its origin 
in the exchange of mesons between the bare nucleons. If one introduces relativity 
into the model, then one has to propose a relativistic Lagrangian describing the 
point-interaction between the nuclei. But this is quite complicated. One of the 
chosen alternatives consists in introducing density dependent energy functionals for 
variational calculations of the Hartree-Fock type. The exchange terms are described 
by ad-hoc ansatze, whose parameters are fitted to experimental data. The classical 
models in this framework are the Gogny and the Skyrme models. These models are 
nonrelativistic. 

Another existing proposition concerns the so-called RMFT {Relativistic mean- 
field theory). The relativistic mean field model is formulated on the basis of two 
approximations: the mean field approximation and the no-sea. On the one hand, 
thanks to the mean-field approximation, the fields for the mesons and the photons 
are treated as classical fields and the nucleons behave as noninteracting particles 
moving in these mean fields. On the other hand, thanks to the no-sea approxima- 
tion, negative energy states belonging to the Dirac sea are not considered and the 
vacuum polarization is neglected. 

The relativistic mean field theory is an effective theory: the Lagrangian of the 
model is an effective Lagrangian with respect to the mean-field and no-sea ap- 
proximations. Since the effective Lagrangian is not derived rigorously from the 
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no-approximation Lagrangian, the parameters of the model must be adjusted on 
experimental data. Therefore, the effects of vacuum polarization as well as the 
correlation effects are not completely neglected but they are taken into account 
implicitly through the adjustment of model parameters. 

During the last years, the relativistic mean-held theory has received a wide at- 
tention due to its successful description of many nuclear phenomena. It has been 
shown that the relativistic mean-held model describes successfully the structure of 
the nucleus and provides a natural explanation for some relativistic effects observed 
experimentally such as the spin-orbit force. This is why the relativistic mean-field 
model can be viewed as the relativistic generalization of some nonrelativistic mod- 
els, such as the Hartree-Fock model with Skyrme or Gogny interaction, where the 
effective forces, which are not appropriate in a relativistic formulation, are replaced 
by average potentials representing independent degrees of freedom. 

For a more detailed description of the RMFT models, see [Sj. The article [I] also 
contains an interesting discussion about these models and their numerical treat- 
ment. 

Although often used in practice, the models of nuclear physics have rarely been 
discussed in the mathematical community. Some nonrelativistic models of nuclear 
physics (of Hartree-Fock type) have been studied by D. Gogny and P.L. Lions in 
an article in 1986 ([2]). To our knowledge, S. Rota Nodari's paper [5] contains the 
first mathematical study of a model from relativistic nuclear physics. 

In the RMFT model studied in this paper, we take into consideration the poten- 
tials created by the mesons a and u>, defining a medium range attractive interaction 
and a short range repulsive interaction respectively. We neglect the meson p de- 
scribing the effects depending on the isospirQ and we omit the photons which are 
related to the electromagnetic interaction (see [5], [3]). This static model which is 
currently known as the a-uj model (|10j. |llj). is given by 

[-ica • V + f3(mc 2 + S) + V] ipj = (mc 2 - fi^tpj , (1.1) 
[-A + m 2 a c 2 ]S = -g 2 a cp s , (1.2) 
[-A + mlc 2 ]V = glcpo, (1.3) 

for j — 1, . . . , A with A the number of nucleons. In these equations the functions 
ipj are the wave functions of the nucleons, c > is the speed of light, m is the mass 

of the nucleons, j3 = y ^ ^ ^ , a& = ( ^ ^ J (k = 1, 2, 3) and the are 

the well known Pauli matrices. The free Dirac operator 

H c = -ica. ■ V + /?mc 2 

acts on 4-spinors, i.e. functions from R 3 to C 4 and it is sclf-adjoint in L 2 (R 3 ,C 4 ), 
with domain iJ 1 (IR 3 , C 4 ) and form-domain H 1 ^ 2 (R 3 , C 4 ). Moreover, ra a and m w are 
the masses of the a and uj meson respectively, and g a and are coupling constants. 



Isospin (contraction of isotopic spin) is a quantum number related to the strong interaction. 
Isospin was introduced by Heisenberg in 1932; he observed that the neutron is almost identical to 
the proton, apart from the fact that it carries no charge. In particular, their masses are close and 
they are indistinguishable under the strong interactions. So, the proton and the neutron appear 
to be two states of the same particle, the nucleon, associated with different isospin projections 
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The densities p s and po are denned by 

A 

Ps =^2^j1pj ■= ^ rj 
3 = 1 3=1 
A 

Po =X^*^'- 

3=1 

In this paper we are interested in the study of the nonrelativistic limit of the 
above system. In order to do so, we are going to pass to the limit as c — > +00 in a 
certain regime for the parameters and we are going to suppose that the "energies" 
of the nucleons are small compared to the speed of light, that is, < pj mc 2 . 



Writing ipj 



Xj 



Ps 



C 2 , the densities p s and po take the form 



A 



I 2 - 1X3 I 2 ) 
1X3 I 2 ) 



3 = 1 
A 

Po =X](N 

3=1 

and the equation ( 1.1 ) becomes 

j -iccr- Vxj + (S+ V)<fj = -pjfj , 

I — iccr ■ V<pj — (2mc 2 + S — V — Pj)Xj = • 



(1.4) 



with er = (a 1, cr 2 , (T3). 

Supposing that the quantity 2mc 2 + S — V — pj is never equal to 0, from the 
system (1.4), we formally obtain 

Xj 



-iccr ■ V(/?,- 



2mc 2 +S-V 



Pj 



(1.5) 



To obtain the nonrelativistic limit of the system (1.4), it would be useful to 
expand the denominator in (1.5) w.r.t. a small parameter when c tends towards 
infinity. In atomic physics, the usual procedure consists in considering that the 
quantity S 2 ^ e 2 Mj is small when c — > +00. In nuclear physics this approximation is 
not acceptable, because the quantity S — V is large. Indeed, for c 
write 



-00, one can 



S = - 



1 



5a 

m„ 



2 r 
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V 



-A 
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Po 
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and, taking into account the physical values of the mesons' masses and of the 
coupling constants (see [5],[S]), we can suppose 

2 / N 2 

+ Ac, 



JJa_ 
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with A > small and 6 > 0. As a consequence, for c — > +oo 

A 

S + V — 29mc 2 \Xj\ 2 - Apo + o(l) , 

A 

S-V = - 29mc 2 \Vj\ 2 + Vo + o(l) , 



and so, from (1.4 1, 



iccr ■ Vxk + ^29mc 2 \xj\ 2 - Xpoj (fk + P-kPk + o(l) = 0, 

t A \ 

iccr ■ Vifk - 2mc 2 - 29mc 2 |</?j| 2 + Xp - p k Xk + o(l) = 0. 

V i =1 J 

Next, if we suppose that = o(l)c-n-oo, we obtain the limiting system 

I A \ 

iccr ■ Vx fe + 29mc 2 ^ IXjl 2 - Xp Pk + P-kPk + o(l) c ^ +oc = 0, 

/ ~, \ 

i<x • Vp k - 2mc 1 - 6} j |<y3j| 2 +o(l) c ^ +00 Xk + o(l) c ^+oo = 0. 
Finally, using the following rescaling 

X 3 = - 2mcV9xj, 
and taking c that goes to +00, we get 

~ \ — ^ 1 - |2 ~ 2ttiX ^ — ^ I ~ |2 ~ , . ~ ~ 

icr • Vx fc + 2^ |Xj I <^fe ij- 2^ I^J I ^ fe + 2m <^ fc = °« 

3=1 J'=l 

- i<r • Vp k l^l 2 j Xfe = 0. 

In this paper we study the above system in the particular case when A = 1, that 
is, the case of a single nucleon. This is of course a very particular case, but it is the 
first step in the study of the multiple particle system which is much more involved 
and that will be the object of a forthcoming paper. The system that we study is 
the following 

( ier ■ Vx+ \x\ 2( P~ a\p\ 2 ip + bp = 0, 

\ -l(T-Vp+(l~\p\ 2 )x = 0, { ""' ] 

<p\ 



with a = ^p, b — 2mp and tp = 

Moreover, we are going to consider a particular ansatz which is classical when 
studying electronic or nucleonic wave functions. It is separable in spherical coordi- 
nates and consists in eigenfunctions of the spin-orbit (angular momentum) operator 
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with corresponding eigenvalue equal to —1 (see [9]). We will look for solutions of 



(1.7) in the particular form 




(1.8) 



where / and g are real valued radial functions. Then the equations for / and g are 

g(f - ag 2 + b) , 

g' = /(i - g 2 ) , 

where we assume /(0) = in order to avoid solutions with singularities at the 



r + 2 f 

r 



origin. For any given value of g(0) there is a local solution of (1.9). Because of the 
fact that we want the corresponding solution if) of ( 1.7) to be square integrable, we 



consider only the solutions of (1.9) which go to (0,0) as r — >• +oo 



About this system our main result is the following 

Theorem 1.1. Given a, b > such that a — 2b > 0, there exists a solution (f,g) € 
C^^+oo^R 2 ) of the system (O) such that /(0) = 0, and th ere exists a constant 
C such that for r > 



with K a i, = min { 



2«- 



2n 



0<-f(r),g(r)<CeM-K a , b r), 



This theorem and its proof have the same flavor as the main results and proofs 
in pQ. In that paper, Cazenave and Vazquez study solutions of the so-called Soler 
model, which consists in a nonlinear Dirac equation. They also consider solutions 
separable in spherical coordinates, with the same angular momentum constraint as 
we do. The main difference between their methods of proof and ours are related 
to the fact that, as we see below, in our case there is a boundedness constraint for 
the initial value of g. This creates additional difficulties and another strategy is 
necessary for the proof. 

The study of the general tridimensional case, without ansatz on the function 
if), involves a completely different method of proof and will be done in a separate 
paper. 



This paper is organized as follows. In Section 2.1 we discuss the possible values of 



the parameters a and b for finite energy solutions of (1.9) to exist. Moreover, some 
preliminary results about the system ( |1.9| are proved there. Section [3] contains 
the proof of our main result, Theorem |1.1| Finally, Section [4] is devoted to the 
presentation of some numerical computations which show the kind of physically 
interesting qualitative properties enjoyed by the solutions obtained in Section [3] 
Those properties will be discussed and compared with existing experimental data. 



2. Conditions on the parameters and auxiliary results 



Proposition 2.1. There is no nontrivial solution of (1.9) such that 
(/(r),ff(r)) 



(0,0) 



as 



(2.1) 



unless a — 2b > 0. Moreover, for all the solutions of (1.9) satisfying (2.1), we have 
g 2 (r) < 1 in [0,+oo). So, in particular, <?(0) must be chosen such that g(0) 2 < 1. 
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To prove the above result, let us first consider the related conservative system 

r/'^-V + O), (2.2) 

This system is the Hamiltonian system associated with the energy 

H(f,g)= 1 -f 2 (l-g 2 ) + ^- b -g 2 . (2.3) 

As a consequence, to have a complete description of the dynamical system, it is 
enough to analyze the energy levels of (2.3). We remind that for every constant C, 
the C-level set is the curve defined by Tc — {(/, g) ; H(f 7 g) = C}. 

Lemma 2.2. For any a, b > 0, H has the following properties: 

(1) if a- b > 0, 

(a) VH(f,g) = if and only if (f,g) = (0,0), (f,g) = (o, iy^ 
(f,g) = (±Vo^M) or if (f,g) = (iv^^-l); 

( b ) (f,g) = f^V^l are local minima ' and (f'S) = (f>g) = 

(±-\/a — b, l) and (/,<?) = (±V a — b, — l) are saddle points of the en- 
ergy H; 

(2) ifa-b= 0,' VH(f, g) = 0if and only if (/, g) = (0, 0) or (/, g) = (0, ±1) ; 

(3) ifa-b< 0, 

(a) VH(f, g) = 0if and only if (/, g) = (0, 0) or (/, g) = ^0, ± yffj ; 

(b) (/, g) — (0,0) and (/,<?) = ^0,±y^j a^e saddle points of the energy 
H. 

Since we are interested in solutions converging to (0, 0) as r — > +oo, we consider 
the zero level set of H. The curve T = {(/, g) ; H{f 1 g) = 0} is an algebraic curve 
of degree 4 defined by the equation 

a 5 4 -2(/ 2 + % 2 + 2/ 2 = 0, (2.4) 

and its behavior depends on the parameters a and b. More precisely, the equation 
(2.4) can be written as 

a (9 - fci,o(/)) (5 + Ko(f)) (g - fta,o(/)) (.9 + ^,o(/)) = , 

with fti, (/) := V , h 2fi (f) := y , provided 

that h (f) = {f 2 + b) 2 - 2a/ 2 > 0. We notice that: 

(1) if a - 2b > 0, then 



ho(f) 



(f - \Ja-b+^aJa-2b)^j (f + y/ a - b + yffia - 2b)j 
(f - \J 'a - b - ^faja - 2b)^j (f + \J 'a - b - ^faja - 2b)^j 
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and ho(f) > if and only if 

/ < -\Ja-b+ y/a(a-2bj, or f > \Ja-b+ y/a(a - 2b) 



or 



-\Ja-b - y / a(a-2b) <f< \Ja-b- v / a(a-2b). 

The zero contour line for this case is represented in Figure [T] 
(2) if a - 2b = 0, then 

h (f) = (/ - V^b) 2 (/ + V^b) 2 > o, 



for all /. In this case the equation (2.4) can be written as 

(bg 2 -f)(g 2 -l) = 0, 

and the zero contour line is represented in Figure [2] 
(3) if a — 2b < 0, then ho(f) > for all /. In this case the zero contour line is 
defined for all / and it is represented in Figure [3j 



=0 o 




-2VT -VT V2~ 2V2~ 




Figure 1. r when a - 2b > Figure 2. T when a - 2b = 




Figure 3. T when a - 2b < 



More in general, Tc is the algebraic curve defined by the equation 
= ag 4 ~2{f + b)g 2 + 2/ 2 - 4C 
= a (g - h ltC (f)) (g + h hC (f)) (.9 - &2,c(/)) (.9 + h 2 ,c(f)) , 



(2.5) 
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with 



'.r-+b+v(f 2 - 


\-b) 2 


- MP 


-2C) 


a 




f by 


- MP 


-2(7) 



and provided that h c (f) := (f 2 + b) 2 - 2a{p - 2C) > 0. By a straightforward 
analysis of the functions h^Ct ^2,c an d using Lemma |2.2[ we obtain the phase 
portrait of the dynamical system ( |2.2[ ) for different values of the parameters a and 
b. In the pictures Figure |2J Figure [5] Figure |6j Figure [7] and Figure [8j the 
grey part represents the region where the energy is negative. 
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Now, we are ready to study the non conservative system (1.9). We begin with 
the following lemma. 

Lemma 2.3. Let x E R. For any a,b > 0, there is t > and (f X i9x) S 
C 1 ([0, t],R 2 ) unique solution of (1.9) satisfying f x (0) = 0, g x (0) = %■ Moreover, 
{fx,9x) can be extended on a maximal interval [0, R x ) with either R x = +oo or 
R x < +oo and limr^fl^ \ f x \ + \g x \ = +oo. Furthermore, (f x ,g x ) depends continu- 
ously on x, uniformly on [0, R] for any R < R x . 

Proof. As in [1], it is enough to write 



s 2 g(s)(f 2 (s)~ag 2 (s) + b) ds , 
r f(s)(l~g 2 (s))ds, 



m 

9(r) 



and note that the right hand side of ( |1.9[ ) is a Lipschitz continuous function of 
(/,(/). The lemma follows from a classical contraction mapping argument. □ 

Remark 1. For any a, b > 0, (f(r),g(r)) = (0,0) is the unique continuous solution 
of (1.9) satisfying /(0) = 0, g(0) = and defined on [0, +oo). 

Remark 2. For a — b > 0, 

1 



9(r)=l 



— yja — b Coth ( Va — b 



r > 
r = 



(2.6) 
(2.7) 



is the unique continuous solution of 1 1.9) satisfying /(0) = 0, g(0) = 1 and defined 
on [0, +oo). Moreover, (—f(r),—g(r)) is the unique continuous solution of (1.9) 
satisfying /(0) = 0, g(0) = —1 and defined on [0, +oo). 
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Lemma 2.4. Let ieR and a, b > 0. Then for any r 6 [0, R x ) we have 

±H(f x (r),g x (r)) = ~/ 2 (r)(l - 5 2 (r)). (2.8) 



Proof. The equation (2.8) is derived by a straightforward calculation 

±H(f x (r),g x (r)) =f x (r){l - gl{r))f x {r) + (-/J(r) + ag x (r) - b)g x {r)g' x {r) 
= -- r f 2 x {r){l-gl(r)). 

□ 

Lemma 2.5. Let a, b > such that a — 2b < and to (f x ,g x ) be the solution of 
(1.9) satisfying f x (0) = 0, g x {0) = x. If x ^ there is no solution that satisfies 
lim (/ a (r),<fc(r)) = (0,0). 

r— >-+oo 

Proof. Let (f x (r), g x (r)) be a solution of ( |1.9| ) satisfying f x (0) — 0, # X (Q) = x an d 
defined on [0,+oo). If lim (f x (r), g x (r)) = (0,0), then \fs > 0, 3S e > such that 

r— >+oo 

\g x (r)\ < e for r > S ej In particular, \g x {r)\ < 1 and the energy H(f x (r), g x (r)) 
is not increasing on (5, +00) for some 5 > (see Lemma 2.4). This implies that, 
Vr > 5, 

ff(/x(r),fc(r))> lim H(f x (r), g x {r)) = 0. 

r— ^+00 

We define by £7 the set 

« = {(/*(r),«fc(r)) G R 2 ; H(f x (r),g x (r)) > and | fe (r)| < l} 
- {(/*(r),0 x (r)) G R 2 ; |<fc(r)| < ft a ,o(/*(r))} , (2.9) 

with 



/i2,o(/x(r)) 



7 2 W + & - V(/* 2 (r) + *0 2 -2a/ x 2 (r) 



Then, (f x (r), g x (r)) € f2 for all r > 5. Moreover, we remark that, if a — 26 < 0, 
then h 2fi (f x {r)) < 1 for all f x (r) € R. 

Suppose now that a; 2 < — . Since H(0,x) < and thanks to the continuity of 
H(f x (r), g x (r)), there exists f < 5 such that H(f x (f), g x (r)) = 0; more precisely let 

f = m{{r € (0,S]\(f x (r),g x (r)) € dn}. 

Since \g x (r)\ = h 2 ,o(f x {r)) < 1, the energy H(f x (r), g x (r)) is not increasing on a 
neighborhood off. As, for all e > 0, H(f x (f - e), g x (f - e)) < = H(f x (f), g x (f)), 
we obtain a contradiction. 

Now, suppose that x > \J^;] in this case H(0,x) > and g x (Q) > 1. Since 
lim (f x (r), g x {r)) — (0,0) and g x is a continuous function, there exists r < S 

r— >-\-oo 

such that g x (f) = 1; more precisely let f = inf{r 6 (0, +oo)|g :c (r) = 1}. This 
implies that <? x (r) > 1 for all r € [0, f), and, as a consequence, H{f x (r),g x {r)) is 
not decreasing on [0, f). Then 

ff(0,a:) <#(/*(?),!) = ^^<0. 
Hence, we obtain a contradiction and g x (r) > 1 for all r 6 [0, +oo). 
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Finally, with the same arguments, we prove that if ac < — y then g x (r) < — 1 
for all r € [0,+oo). ' □ 

Lemma 2.6. Let a, b > swc/i £/ia£ a — b > and let (f x ,g x ) be the solution of 



1.9\ ) satisfying f x (0) = 0, g x (0) = x. If x 2 > 1, then g 2 (r) > 1 for all r € [0,R X ). 
Proof. To prove this lemma, we use the monotonicity properties of the function 

a 4 _ b r 
4 . 2" 

and increasing on I +oo 
Let g x (0) = x > 1 and suppose, by contradiction, that there exists rg such that 



F(x) — |x 4 — \x 2 . In particular, we use the fact that F(x) is decreasing on 



9x(fo) = 1 and g x (r) > 1 for all r G [0,r ). As a consequence of Lemma 2.4 the 
energy H(f x (r),g x (r)) is not decreasing on [0,r ), which means 

H(0,x) < H{f x (r ),g x (r )) , or equivalent^, F(x) < F(l). 



Since y ^ < 1, the above inequality contradicts the monotonicity of i 7 ". As a 
conclusion, ^(r) > 1 for all r G [0, R x ). In the same way, using the fact that 
energy is not decreasing if g x (r) < —1, we prove that if g x (0) — x < —1, then 
2 x (r)<-l for all re [0,^). □ 

Lemma 2.7. Let a, b > such that a — 2b and let (f x ,g x ) be the solution of 



1.9) satisfying f x (0) ~ 7 ^(O) = x. If x ^ there is no solution that satisfies 



hm (f x (r),g x (r)) = (0,0). 



Proof. Let (f x (r), g x (r)) be a solution of ( |L9| ) satisfying f x (0) = 0, ^(0) = x 
and defined on [0,+oo). Suppose first that x 2 > 1. Then, thanks to Lemma 



2.6 g x (r) > 1 for all r 6 [0, +oo). Next, suppose that x 2 < 1. As before, if 
lim (f x (r),g x (r)) = (0,0), then Ve > 0, 3S e > such that \g x (r)\ < e for r > 

r— >-+oo 

S £ . In particular, | (r) | < 1 and the energy H(f x (r), g x (r)) is not increasing on 
(<5, +oo), for some 5 > 0. This implies that, Vr > 5, 

H(f x (r),g x (r))> lim H(f x (r), g x {r)) = 0. 

Since i?(0, x) < and thanks to the continuity of H(f x (r), g x (r)), there exists f < 5 
such that H(f x (f), g x {f)) — 0; more precisely let 

r = mi{re(0,5]\H(Ur),g x (r))=0}. 

By definition of f, we have that g 2 (r) < 1 for all r € [0,r); this implies that the 
energy is not increasing on [0, f). Then 

= H(f x (r),g x (r))<H(0,x), 

a contradiction. Finally, thanks to Remark [5J if a; 2 = 1, then g 2 (r) = 1 for all 
r e [0, +oo). □ 

Lemma 2.8. Let a, b > such that a — 2b > and let (f x ,g x ) be the solution of 



1.9) satisfying f x (0) = 0, g x (0) = x. If x 2 < 1, then gl(r) < 1 for all r G [0,R X ). 



Proof. To prove this lemma, we use the monotonicity of the function F(x) — |x 4 — 

In particular, since W 2 ^ < 1, 



\x 2 and the fact that F(x) < in 
F(x) < F(l) for all x such that x 2 < 1 
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Let g x (0) = x such that x 2 < 1 and suppose, by contradiction, that there exists 
r such that g 2 (ro) = 1 and g 2 (r) < 1 for all r G [0,r ). As a consequence of 
Lemma 2.4 the energy H(f x (r), g x (r)) is not increasing on [0, ro), that means 

H(0,x) > H{f x (r ),g x (r )) , or equivalently, F(x) > F(l). 

The above inequality contradicts the properties of F. As a conclusion, g 2 (r) < 1 
for all r G [0,R X ). □ 

Proposition 2.9. Let a, b > such that a — 2b > and fe£ g x ) be the solution 
of (1.9) satisfying f x (0) = 0, g x (0) = x and such that lim (f x (r), g x {r)) = (0,0). 

T— 5- + 00 

Then x 2 < 1 and g 2 {r) < 1 for all r € [0, +oo). 



Proof. Let (f x ,g x ) be the solution of ( |1.9[ ) satisfying / x (0) = 0, g x {0) — x and such 
that lim (f x (r),g x (r)) = (0,0); then Ve > 0, 3S £ > such that \g x (r)\ < e for 

r > 6 S . Hence, as a consequ ence of Lemma 2.6 and Remark[2j we have x 2 < 1. 



Finally, thanks to Lemma 2.8 we can conclude that g 2 (r) < 1 for all r G [0, +oo). 

□ 



Proof of Proposition \2. 1\ It follows readily from the above series of lemmata. □ 

3. Existence of solutions 

In this section, we prove Theorem |1.1| that is, we show that whenever a > 2b, 
there exists at least one nontrivial solution of the system 

(31) 

s' = /(l-9 2 ), 

which goes to (0, 0) as r — s- +oo and such that /(0) = 0. Moreover / < and 
< g < 1 (or equivalently, / > and — 1 < g < 0). 

In what follows, let then a, b > be such that a — 26 > and let (/a,, g x ) be the 
solution of (3.1 1 satisfying ./^(O) = 0, (fc(0) = x. 

We remind that (3.1) is associated to the conservative system 

jf'=9(f 2 -ag 2 + b), 

and we define A = {(/o,5o) G R 2 | 2/ 2 - - (a - 2b) < 0, g 2 < 1}, the set of 
admissible initial conditions. We recall that the energy of the system is given by 

H{f,g) = \f{l-g 2 ) + -/-\g 2 . 



a-2b 
4 ' 



Remark. If g 2 < I, (f ,g ) € A if and only if H(f ,g ) < H(0, 1) = 

Lemma 3.1. Let (/o,<7o) G A and let (/,(?) be the solution of the conservative 
system (S.Ofy with initial data (/o,5o)- Then (f(r),g(r)) G A for all r G [0,+oo). 
In particular, (f(r), g(r)) is bounded for all r G [0, +oo). 

Proof. First of all, let (fa, go) G A such that g 2 < 1. Since in the case of autonomous 
systems two different trajectories cannot intersect, if g 2 , < 1 then, g 2 (r) < 1 for all 
r G [0,+oo). Moreover, H(f,g) = H(f ,g ) < 2=^. This implies 

2f 2 (r)-ag 2 (r)-(a-2b) <0, 
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for all r G (0, +00). 

Now suppose that go = 1 (go = — L respectively). In this case, ^(r) = 1 
(g>(r) = —1, respectively) for all r G (0, +00) and /(r) solves 

//'(r) = / 2 (r)-a + & ) 

1 /(0) = /o • ■ (3 - 3) 

Hence, if /o = ±V a — b, then /(r) = ±Va — & for all r G [0, +00). Finally, if 
fo 7^ ±%/a — 6, we solve the Cauchy problem (3.3) and we obtain 

M = -2^1 (l - - 1)"' - v^. 

Note that f(r) is bounded for all r G (0, +00). Indeed, if (fo, 1) G .4 and / 7^ 
±Vct — 6, then — ya — < fo < ^/a — b and /'(r) < for all r G (0, +00). Hence, 
/o > f(r) > -V^b for all r G (0, +00). □ 



Lemma 3.2. Let (/o,5o) G -4- (/,<?) 6e i/ie solution of (3.2) with initial data 
(/o,<7o)- Let (f„,gn) € -4 Pn &e sitc/i that 

lim p„ = +00 and lim (/°, #°) = (/o, So)- 

n— f+oo n— f+00 

(fnidn) be a solution of 

2 

/n H — /n = ffnC/n ~ affn + &) . 

Pn + r 

5^ = fn(i- gl), 

such that /n(0) = /„, g n (0) — <?„) an d let [0j T n) be the maximal existence interval 
of(f n ,g n )- Then 

(1) lim T n = +00, 

n— f+oo 

(2) (fni9n) converges to (/,<?) uniformly on bounded intervals. 

Proof. For the proof of this lemma, we proceed exactly as in the proof of Lemma 
2.5 of pQ; the only difference is that we do not work in all R 2 but in A. For the 
reader's convenience, we rewrite the proof in our case. 

We have to prove that for any R < +00, lim r n > R and (f n ,gn) converges 

to (/, g) uniformly on [0, R]. Then consider R < +00. Thanks to Lemma 3.1, if 
(/oj go) € then (J (f(r),g(r)) is a bounded set. Hence, applying a contraction 

re[0,jR] 

mapping argument, we obtain that there exists S > and t > satisfying the 
following property: for any p > 1, tq € [0, i?] and (wo,^o) G -4. such that |u>o — 
/( r o)l + l z o — gi^o)] < there exists a solution (w,z) G C 1 ([0, r],R 2 ) of 

2 

u/ H w = z(w 2 — az 2 + b) . 

P + r (3.4) 

z' = w(l-z 2 ), 

with iu(0) = wo, z(0) = 20 an d SU P |w(r)|+|z(r)| < 2M, where M = sup |/(r)| + 

re[0,r] r-e[0,r] 

W)\. 
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Next, let K the Lipschitz constant of the second term of the system (3.4) on the 
ball of radius 2M. We get, for all r £ [0,r], 

\w(r)- f(r + r )| + \z(r) - g(r + r„)| < ho - /(r )| + No - ff(r )| 
2 



|w(s)|ds+ / if(|t«(fl)-/(* + ro)| + |2!(«)-ff(* + ro)|)dfl 
P + s Jo 

4Mr 

<K " /Ml + No - SMI + 

P 

+ f K(\w(s) - f(s + r )| + \z(s) - g(s + r )|) da. 
Jo 

Thus, thanks to Gronwall's lemma, we have 

4Mt^ 

for all r € [0, t]. Finally, we apply this argument with p = p n , r$ = 0, wq = /„, zo = 
<7°. We obtain lim r„ > r and (f n ,gn) converges uniformly to (/,#) on (0, r). To 

conclude, we iterate the above procedure N times with Nt < R < (N + l)r. □ 

Lemma 3.3. Let a, b > such that a — 26 > and Ze< (f x ,g x ) be the solution of 
(3.1 ) satisfying f x (0) = 0, g x (0) = X. If x 2 < 1. then gl(r) < 1 and fl(r) < a - b 
for all r £ [0, R x ) and R x = +oo. Moreover, (f x (r), g x (?)) € A, for all r £ [0, +oo). 



\w(r)-f(r + r )\ + \z(r)-g(r + r )\ < \w - f(r )\ + \z - g(r )\ + 

P 



e K ^ 



2.4 



Proof. By Lemma 2.8 we have g x (r) < 1 on [0,R X ). Then, applying Lemma 
we obtain that the energy is non-increasing. Thus, 

H(f x (r),g x (r))<H(0,x) < , Vr e [0, R x ) , 

and by the remark following the definition of the set A, (f x (r),g x (r)) £ A and 
fx( r ) < a — b, for all r £ [0, R x ). In particular, R x = +oo. □ 

Lemma 3.4. Let a-2b>0,0<x<l and f x (r) < 0, g x (r) > for r £ (0, +oo). 
Then there exists a positive constant C such that 

< -f x {r),g x (r) < Cexp(-K a , b r) , (3.5) 

with K a<b = min{|,^}. 

Remark. The ideas behind the proof of this lemma are inspired by the work of 
Cazenave and Vazquez pQ. 

Proof. Applying Lemma 3.3 we have g x (r) < 1. Thus g x is decreasing on (0, +oo) 
and there exists 6 > such that lim g x (r) = S < 1. 

r— v+oo 

Suppose for the moment S = 0. Then lim g x (r) = and g x (r) < i- for r large 

r > • oc 

enough. Considering (3.1) we have 

f x >--f x + ^g x , g x : , (( /,. 

Hence, for r large enough, 

(flx-/ B ) / +^«,6(ffx-/*)<0 ) 

with -fT^f, = mm{|, 2 % }■ Integrating this inequality, we obtain 
~fx{r) + g x (r) < Cexp(-iT 0)6 r). 
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Therefore it only remains to prove that 6 = 0. By contradiction, suppose S > 0. 
Let {r n } n be a sequence such that lim r n = +00 and lim f x (r n ) — Let 



[u, v) be the solution of ( |3.2| with initial data (k,5). It follows from Lemma 3.2 
that (f x (r + •)) converges uniformly to (it, v) on bounded intervals. Since 

lim g x (r n + r) = 5 for any r > 0, we have u(r) = 5 < 1 for any r > 0. Hence, 

n— >-\-oc 

from the second equation of (3.2) and applying Lemma 2.2 we have k = and 
S = \ -. As a conclusion, this means 



Urn (f x (r),g x (r))= 0, ■ 



r— >-+oo 



Now, we put (fc 



to. First of all, we remark that 



lim (f x (r)Mr)) = (0,0) 

r— > +00 



and w > on [0, +00). Furthermore, the equations for (/ x , w) are given by 

2 

/■; + ^ 




which implies that as r — > +00, 



w 



-2bw + o(\w\ + \f x \) 
b 



U 1 



(3.6) 



+ o(H + l/*l). 



Then, 



1 /a 
2 



(«> + /*) + 



a - b 2 
r 



2n 



2o 



+ 26 ttf + odwl + l/xl). 



As a consequence, for r large enough, we obtain 

b 



(w + f x )'- \ '" 



(w + f x )<0. 



Therefore, exp (- 



[w + f x ) is non-increasing for r large, and since 



lim exp 

r— »+co 



ft — 6 

2a 



r)(w + f x )=0, 



we should have w > —f x > for r large. This implies that o(\w\ + \f x \) can be 
replaced by o(|u?|). To conclude, we put f x = —f x and we rem ark t hat f x > 
on (0,+oo) and w > f x for r large. From the first equation of (3.6) we obtain, 
f ) fx > 6/ x , which implies /^(r") > Ce br for some C > 0, 



'36 



for r large, > , ., 
contradicting lim f x (r) = 0. Thus the lemma is proved. 



□ 
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We are now ready to prove the main result of this paper. The proof of this 
theorem is given in several steps. Some parts of this proof are inspired by the work 
of Cazenave and Vazquez [I]. For our problem, we have to introduce new ideas to 
deal with the fact that the initial condition for g have to be bounded. 

We define the set / as follows: 

/ = lx > y| | 3r x e (0, R X )J X < and g x > on (0, r x ), f x (r x ) = j . (3.7) 

In the first step, we show that / is a non empty open set. In the second step, we 
prove that / is a bounded set and that sup/ < 1. Finally, in the third step, we 
conclude taking x = sup J as initial condition. 

Since we take x — sup / as initial condition for g, it follows from Proposition |2.1| 
that sup / must be strictly smaller than 1. Actually this is the main difficulty to 
deal with and represent the principal difference between our proof and Cazenave- 
Vazquez's one. 

Step 1. / is a non empty open set. 

To prove that / is a non empty set, we show that I \ \ — ) C /. Indeed, let 



x G I — I C I ). From Lemma 



3.3 



we have g x (r) < 1, f x (r) < a — b 



and R x = +oo. Hence the energy is non-increasing and 

H{f x {r),g x {r)) < H(0,x) < H (o, =0, 

for r > 0. On the other hand, H(y, 0) = ^y 2 > for y e R. As a consequence, g x 
cannot vanish in a finite time. Therefore g x (r) > for all r £ (0, +oo). 



Next, using the first equation of (3.1 ), we obtain f x (0) = x(b—ax 2 )/3 < 0. Hence, 
fx{r) < for r > small enough. Thus it remains to prove that f x (r) vanishes for 
some r > 0. Suppose by contradiction that f x (r) < for all r e (0, +oo). Then, 
we can apply Lemma 



3.4 



which asserts that lim (f x (r), g x (r)) = (0,0). This 

r— »+co 

implies lim H(f x (r), g~x[r)) — 0, which is a contradiction since H(f x (r), g x (r)) < 

T— 5- + 00 

Z/(0, x) < 0. Therefore f x has to vanish and x € /. 

The open character of / is a consequence of the continuous dependence of (f x , g x ) 



Step 2. sup / < I. 

First of all, we show that sup / < 1 and we state some properties of solutions of 



(3.1 ) with initial condition x G /. 



Lemma 3.5. Let x £ I. Then < g x {r x ) < y - and H(f x (r x ), g x (r x )) < 



Proof. Suppose by contradiction g x (r x ) > y Using the first equation of (3.1), 
we get 

f' x (r x ) = g x {r x ){b - ag 2 x (r x )) < 0. 
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Hence, f x is decreasing in a neighborhood of r x ; this implies that for e > suffi- 
ciently small f x (r x — s) > f x (r x ) — 0, which contradicts the definition of r x . Finally, 

if < g x {r x ) < Jj, then H(f x (r x ),g x (r x )) = %gi{r x ) - \g%{r x ) < 0. □ 



As a consequence of Lemma 



3.5 



I C I \ -, f . Indeed, if x > 1 then, thanks to 



Lemma 2.6 and Remark[2j g x (f) > I for all r G [0, R x ). Hence sup/ < I. 
Lemma 3.6. Let x e I. Then < -f x (r) < ^f\g x {r), for all r £ (0,r x ). 

Proof. Let h x (r) = ^/%g x (r) + f x (r) for < r < r x . We have h x (0) = x > 
and h x (r x ) — 9x(y x) > 0. Suppose by contradiction that there exists p € (0, r^) 



such that h x (p) < 0. Then, there exist pi,p2 € (0,r x ) such that 
pi < p, hx(pi) = and ft^(pi) < , 
92 > P, h x (p 2 ) = and h' x (p 2 ) > 0. 

Next, we remark that 

f^) = (^ + |^)-/»-^« r) 



Hence, 



(/|(r-)-a^( r ) + b )^( r ). 



a - 2b 2 fa 



2 p! V 2 
and, since f x (pi) < 0, we obtain 

2\/2a 



fx(pi 



Pi 



> 



(a - 2b) ' 

In the same way, h' x {p2) > implies 

2V2a 



P2 < 



(a - 2b) ' 
As a consequence, 



2V2a 

a contradiction. □ 

Let us next prove that the supremum of / cannot be equal to 1. For this, we use 
a contradiction argument. We start with some estimates for solutions of (3.1 1 with 
initial data x n £ I supposing that lim x n = 1. 

n— >-+oo 

Lemma 3.7. Let {x n } n <Z I a sequence of initial conditions such that lim x n = 

n— v+oo 

1, and define j Xn = min{r > : f' x (r) = 0}. Then there exists a constant c > 
such that 

inf j x = c. (3.8) 

raGN 
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Remark 3. Since f' x (0) < 0, j Xn is such that 

f' Xn (r) < on [0,7, J and f' Xn (j Xn ) = 0. 
As a consequence, f Xn ( , y Xn ) < 0. 

Proof. For convenience we denote f Xn , g Xn and j Xn by /„, g n and 7„. 

Since x n — > 1 as n — > +oo, let us assume that i n is such that < x 2 n < 1. 
Then there exists S Xn e (0, +oo) such that 

9n(S Xn ) = x 2 n and < #„(r) < x n Vr e [0, S^J . 
Let denote 5^ by S n . Hence, 

x n - xl = g n (0) - g n (S n ) = - [ /„(r)(l - gl(r)) dr = [ ( k n (s)s 2 ds^ , 

Jo Jo Jo r 

with 

K(s) = 9n (s)(l g 2 n (s))(f n (s) + agl(s) b). 

Next, we remark 

x n (l - x n )m afi < k n (s) < x n (l - x n )M afi , (3.9) 

where m a ^ b and M a ^ are strictly positive constants that depend only on the pa- 
rameters a and b. As a consequence, 

x n (l-x n ) rn ^ b -S 2 < [ f k n (s)s 2 ds ^ < x n (l - x n ) ^^Sl , 
o Jo Jo r o 

which implies 

~M n h r 1/2 . „ . rm 8 , 6 i-i/2 



< 5 n < 

Hence, there exists a constant c > such that inf„ e ^ S n = c. 

Now, suppose by contradiction that inf„ elN | 7„ = 0, then there exists a subse- 
quence of {j n }n such that 7„ — » and, for n sufHciently large, 7„ < S n . 

n 

First of all, we remark that f n {ln) satisfies 
2 

— fnhn) = 9nhn)(fn(ln) ~ C^n^n) + b) , (3.10) 
7n 

and it has to be strictly negative. Furthermore, since x 2 n < g n {ln) < x n , for n large 
enough we can write 



As a consequence, 



lim fniln) = lim \—-\\+ 9l{ln){agl{"f n ) - b) ) 1 = 0. 

™^+°° n^+oo \^7„ y 7^ J g n {in) 

By the definition of j n the function /„ is decreasing on [0,7„). Then, 

> f n (r) > / n (7„) Vre(0, 7 „), 

and, passing to the limit, we obtain that for all e > there exists n (e) G N such 
that 

l/n(r)| <e, 
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for all n > no(e) and for all r € [0, j n ]. Moreover, using the second equation of 
the system ( |3.1| ), we have that g n is decreasing on [0, 7„). Thanks to the fact that 
x n < 9n{ln) < in, we get for all e > there exists mo(e) € N such that 

l-e<0n(r) <l+e, 

for all n > m (e) and for all r £ [0,7„]. Hence, for all e > and for n sufficiently 
large, we have the following estimates for all r £ (0, j n ] : 

(1 - e)(b - a(l + ef) < f n {r) + -J n (r) < (1 + e)(b - o(l - ef + e 2 ) , 
which implies 

(l- £ )(6- a (l + £ )V 2 ^ (l + £ )(&-a(l- £ ) 2 + £ 2 )r 3 

and, as a consequence, 

(l- £ )(b-q(l + £ ) 2 ) < / n ( 7n ) < (l + £ )(fe-a(l- £ ) 2 + £ 2 ) ^ 
3 ~ 7n ~~ 3 



for all e > and for 7i sufficiently large. Using (3.11) and the first equation of the 
system ( |3.1[ ), if e is chosen sufficiently small we obtain 

2 

/n(7n) = fn{ln) + 9n{ln){fl{ln) ~ ag n (j n ) + b) 

In 

2(l- £ )(6-a(l + £ ) 2 ) , ,,, , ,„ 2 , b-a „ 
< — + (1 + e)(b - a(l - £ ) 2 + £ 2 ) < — < , 

which contradicts the definition of 7„. So, there exists a constant c > such that 



inf 7„ = c. 

new 



□ 



Lemma 3.8. Let {x n } n d I a sequence of initial conditions such that lim x n = 

n— >+oo 

1, and define j Xn — min{r > : f' x (r) = 0}. Then, up to a subsequence, 

lim 73; = +oo. (3-12) 

n— >+oo 



Proof. Thanks to Lemma 3.7 we know that the sequence {"f Xn } n is bounded from 
below by a constant c > 0. 

Now, suppose by contradiction that sup„ eN 7 Xll — C < +oo. Using the conti- 
nuity of the flow and the fact that the sequence {"f Xn } n is bounded from above by 
some constant C, we obtain 

V £ >0,3n £ >0,Vne N,n>n e => \gi{j x J - g Xn (7x n )\ <e, 
which means 

lim 5 , x„(7x„) = l, (3.13) 

n— J-+oo 

since = 1 for all r > 0. With the same arguments, we get 



lim 

n— >-+oo 



(~ - v^&Coth (Va = &7*„) " /x„(7*„)J - (3.14) 
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where ^ \fa — 6Coth (%/a — b-f Xn ) is the expression of /i(7x n )- Next, we 

remark that f Xn ("f Xn ) satisfies 

■ 2 /x n (7x n ) = ffx n (7xj(/x„(7x„) - a<?x„(7x n ) + b) , (3.15) 
and it has to be strictly negative. Hence, for n large enough, we can write 



/x n (7xj= (^--4/4-+^ n (7*»)(o52 B (7x n )-6)J -4r^- ( 3 - 16 ) 
7z„ y 7a; 7l / 9x n \1x n ) 



So, using (3.14), we have 

lim ^J-^- + (a^b)- y/a - bCoth (y/a - bj Xn ^ = 0. (3.17) 

Let /i(r) be defined by 

h(r) := + (a-b) - Va- 6Coth (-y/a-fcr) . (3.18) 
/i is a smooth function such that lim h(r) = — lim ?i(r) and /i(r) > for 

r— 7-0 r— v+oo 

r e (0, +oo). Indeed, 



+ (a - 6) > Va-6Coth(Va-6r) 

if and only if 



1 



r 



+ (a - 6) > (a - b) Coth 2 (Va - br) , 



or equivalently, if and only if 

Sinh(-\/ a — b r) > \J a — br . 

Let us analyze the function \{x) = Sinh(x) — x. The function l(x) has the following 
properties: 

(1) lim l{x) = 0; 

(2) V(x) = Cosh{x) - 1 > on (0, +oo). 

As a conclusion, / > on (0, +oo) and h(r) > for r E (0, +oo). Furthermore, 
thanks to the fact that the sequence {7x n }« is bounded from above, up to a subse- 
quence, 

lim j Xn = 7 



with 0<c<7<C< +oo; hence we obtain 



lim /i(7x„) = /i(7) > 0, 



which contradicts the limit in (3.17). As a consequence, sup neN 7 Xji = +oo and 
there exists a subsequence of {"f Xrl }n such that 



lim j Xn = +oo. 



□ 
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To summarize, let {i„}„ C / be a sequence of initial conditions such that 



lim x n = 1, 7 X = min{r > : f' x (r) = 0} and S Xn = min{r > : g Xn (r) = x 2 n } 

n— f +oo 

Then, up to a subsequence, lim j Xn = +00, 



as in the proof of Lemma 



inf S Xn > and sup S Xn < 



3.7 



00. Moreover j Xn < r Xn . 
Let us now end Step 2 by proving the following lemma. 
Lemma 3.9. Let x = sup/. Then x < 1. 



Proof. For convenience we denote <7 Xn , 7 Xn and by /„, 7„ and 5 n . We 
already proved that sup / < 1. Then suppose, by contradiction, that there exists a 
sequence {x n } n C / such that lim x n = 1. 

We claim then that in this case there exists a sequence {t n } n such that S n < 
t„ < 7« for all n £ N and lim \f n {t n )\ + \g' n (t n )\ = 0. 

n— >-\-oo 

Indeed, suppose by contradiction that, up to a subsequence, 
lim inf \f' n (r)\ + \g' n (r)\>0- 

n->+oo 7'6[S„,7„] 

We define 7V x (i?i,i? 2 ) as the number of roots of g x in {r G (RijRz) ; g^(r) + 
/J(r) 7^ 0} and use the notation N X (R) = N X (Q,R). We remark that the number 
of roots N x (Ri, R 2 ) is linked to the winding number of the trajectory around (0, 0) 
in the following way: 

W|fl2) = WzW i (319) 

where 6 x (r) — — arctan ^ g*^) ) ■ Then, we compute the derivative of 8 x (r) and we 
get 

We note that 6*^(0) = > °- Hence ^( r ) 

is increasing in a neighborhood of 

0. Let denote N Xn and 9 Xn by N n and 8 n . 

First of all, we remind that for all n £ N the following properties hold for all 

r £ [S n ,j n ]: 

• f'n(r) < 0; 

• — \fa — b < fni'Jn) < fn( r ) < fn{S n ) < S < for some S > that does 
not depend on n; 

• g'Jr) < 0; 

' \pk < 9n(ln) < 9n(r) < 9n(Sn) = X 2 n < 1. 

Note that the first inequality of the last property is obtained using the first equation 
of (3.1 1 and the definition of 7„. Moreover, the inequality f n (S n ) < — 5 < will be 
proved below. 

Using these properties and the definition of N n (S n , j n ), we obtain 

N n (S n ,i n ) = (3.20) 

for all ri £ N. 

On the other hand, as we prove below, we have 

lim iV n (5 n ,7„) = +oo. (3.21) 

n— J-+oo 
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Indeed, 

1 f 7 " 

lim inf N n <S n , 7«) = lim inf B (r) dr 



> — liminf inf 9' n ( r ) ( um 7n + hm sup SV. 

7T n->+oo re[S„,j n ] ^n-^+oo n^+oo 

Hence, if liminf inf 6^(V) > 0, we obtain 

n^f+oo rG[S„,7„] 

liminf A^„(7 n ) = +oo 
since lim 7„ = +oo. Therefore, it remains to prove 

liminf inf 6' (r) > 0. (3.22) 

rn-+oo rG[S„,7„] 

We remark that, for all r £ [5 n ,7„], 

/»W„(r) - /;(r) fln (r) |/„(r)|K(r)| + |/4(r)|| 9n (r)| 



#(0 + <£(r) /S(r) + «£(r) 

min < <5, 



which implies 



a — o + 1 



inf ^(r)>C inf \g' n (r)\ + \f' n (r)\ 

i-e[S„,7Ti] re[S n ,7„] 



min{5,vT) 

with C = — a _ b+ i > °- Hence, 

liminf inf 0' n (r) > C lim inf \g' n {r)\ + |J»| > 0. 

n^+oor6[S„,7„] ri->-+oo r g [S„ ,7„] 



As a conclusion, if up to a subsequence lim inf \fL( r )\ + ls«( r )l > 0, (3.20) 

n->+oo re[S„,7„] 

and (3.21| provide a contradiction. 

Hence, there exists a sequence {i n }„ such that for all n € N, SVi < in < 7n and 

lim |/;(*n)l + 1^(01=0. 

n— > +oo 

With the same arguments used in the proof of Lemma |3.8[ we obtain, up to a 
subsequence, 

lim t n = +oo. 

n— f+oo 

Indeed, the only difference between this case and the previous one lies in the fact 
that f'('Jn) — while f'(t n ) tends only towards as n goes to +oo. Then, the only 
change to be made in order to adapt the arguments used in the proof of Lemma 



3.8 to the present case consists in adding a small o(l) n _ > + 00 to the r.h.s. of both 



( |3.15| and ( |3.16 |. 

Next, we remark that 



limsup/ n (t n ) < limsup/„(S'„) < 0. (3.23) 

n— >-+oo n—>-\-oo 
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Indeed, since S n < t n < j n for all n G N and by the definition of j n , we have 
fniln) < f n (t n ) < fn{S n ) for all 71 G N. It remains to prove that limsup f n (S n ) < 



n— >-+oo 



0. Using the second equation of the system (3.1 ), we have 



Snfn(Sn) < / fn(f) dr 



9'n(r) 
l-flS(r) 



dr 



1, 1 

2 log 7T 



< s < o, 



< 1 uniformly. As a consequence, and using the fact that {5 n } is 



since ( 1+Xn )2 
bounded, we find 

lim sup /„(£„) < 0. 

n— )-+oo 

Therefore, by the equation satisfied by <?„, lim g' (t n ) = implies lim g n (t n ) = 

n— >+oo n— f+oo 

1. Hence, using the first equation of (3.1) and the properties of convergence of 
{/«.(*«)}« an d {^n}n, we obtain, up to a subsequence, lim = a — b. As a 

n— f+oo 

conclusion, lim f n (t n ) = — V a — b since f n (t n ) < 0. 

n— >+oo 

Next, let us set (w„, z„) = (/„(*„ + -))5n(*n + ■))• We nave 



lim ui„(0) = — v 7 a — b , 
n— ^+oo 

lim z n (0) = 1. 

n— >+oo 



Let (it?, 2) be the s oluti on of (3.2 1 with initial data (—y/a— b ,1). Since (— ya — b , 1) 
is a rest point of (2.2), w(r) = — ya — 6 and z(r) = 1 for all r G [0, +00). Hence, 
thanks to the condition a — 26 > 0, we obtain 



(3.24) 



3.2 



and since lim t Tl 

71— >+00 



for all r G [0, +00). On the other hand, applying Lemma '• 

+00, we have that for all r G (0, +00), (w n ,z n ) converges to (w,z) uniformly on 
[0, r]. Therefore for n large enough and for all r G (0, +00), we have 



□ 



fn(t n + r ) + \l ^gnitn + r) < , 



which contradicts Lemma |3.6| This proves the lemma. 
Step 3. Conclusion and proof of Theorem 

Let (/,<?) = (fx, 9x) with x = sup I. Th anks to Lemma 



oil. 



< x < 1. Hence, by Lemma 



3.3 



3.9 



and by the definition 



we have Rx 



-00. Since x > 



have /'(0) < 0, and thus / < and g > in a neighborhood of 0. Moreover, since 
the set I is open, x ^ I and / cannot vanish before g does. 

On the other hand, x G I. Then, g cannot vanish before / because of the 
continuous dependence of (f x , g x ) with respect to x. F inally, / and g cannot vanish 
simultaneously because (0,0) is a rest point of (3.1). Therefore, f(r) < and 
g(r) > for all r G (0, +00). Then, using Lemma 3.4 there exists a constant C 
such that for all r G (0, +00), 

< -f x (r),g x (r) < Cexp(-K a , b r) , 

withX Q)b = mm{£,^ i '}. □ 
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4. Numerical aspects and qualitative properties of ground states 

As it has been observed in computations carried out in physics, the fields of 
mesons a and u) are proportional to the scalar and the vector density. In models 
for finite nuclei these fields approximately assume a "plateau" -like Saxon- Woods 
shape: they vanish outside the nucleus and they are more or less constant inside it. 
Moreover, the intensity of the potential for the antinucleons V — S is much more 
higher than that of the potential for the nucleons V + S (see for instance [5]). 

We have run some numerical calculations for our model, trying to see how the 
values of the parameters in the problem affect the shapes of those fields and the 
intensity of the potentials. We observe that the results depend a lot on the values 
of a and b, which are related to the physical values of the meson masses and of 
the coupling constants. More precisely, we remark that even in this very particular 
case where only one nucleon is taken into account, the Saxon- Woods shape for the 
potentials V and S is perfectly observable and the magnitude of \V — S\ is much 
larger than that of |V + S\, if the parameters a and b are well chosen. 





Figure 9. Plot of a 2 , ^ ~ 1 

& 5 a 



Figure 10. Plot of/ 2 , 



2/j 



In Figures [9] and 10 we plot the shapes of f 2 and g 2 for the values a = 9 and 
6 = 4. 

The "plateau" -like Saxon- Woods shape is perfectly clear for g 2 . For f 2 it is 
not, but taking into account that as c goes to +oo, S ~ —mc 2 g 2 + f 2 /(4m) and 
V ~ mc 2 g 2 — ag 2 1 (2m) + f 2 / (4m), the properties discussed above for V, S, V + S 
and V — S are verified here in a very straightforward way. Indeed, the above 
asymptotics show that V and S behave like a plateau if g 2 does, and the intensity 
of \V - S\ is much higher than that of \V + S\. 

A less clear case is the next one, see figures [TT] and [l2j where the values of a and 
b are respectively of 4 and 1. Here the "plateau" is much less visible and its edge 
are much less sharp. Practically there is no "plateau" in this case. 



1 



S 2 (r) 




Figure 11. Plot of g 2 , ~ ^ 1 




Figure 12. Plot of / 2 , f ^ 1 
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Actually, other computations that we have run show that the "plateau" shape 
is more and more present when 2b /a approaches 1. 
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